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Abstract: Light stops consistent with the Higgs boson mass of ^ 126 GeV are investigated within 
the framework of minimal supergravity. It is shown that models with light stops which are also 
consistent with the thermal relic density constraints require stop coannihilation with the neutralino 
LSP. The analysis shows that the residual set of parameter points with light stops satisfying both 
the Higgs mass and the relic density constraints lie within a series of thin strips in the mo — mi /2 
plane for different values of Ao/mo. Consequently, this region of minimal supergravity parameter 
space makes a number of very precise predictions. It is found that light stops of mass down 
to 400 GeV or lower can exist consistent with all constraints. A signal analysis for this class 
of models at LHC RUN-H is carried out and the dominant signals for their detection identified. 
Also computed is the minimum integrated luminosity for 5a discovery of the models analyzed. If 
supersymmetry is realized in this manner, the stop masses can be as low as 400 GeV or lower, and 
the mass gap between the lightest neutralino and lightest stop will be approximately 30-40 GeV. 
We have optimized the ATLAS signal regions specifically for stop searches in the parameter space 
and find that a stop with mass ^ 375 GeV can be discovered with as little as 60 fb ^ of integrated 
luminosity at RUN-H of the LHC; the integrated luminosity needed for discovery could be further 
reduced with more efficient signature analyses. The direct detection of dark matter in this class of 
models is also discussed. It is found that dark matter cross sections lie close to, but above, coherent 
neutrino scattering and would require multi-ton detectors such as LZ to see a signal of dark matter 
for this class of models. 


1 Introduction 

The discovery of the Higgs boson in 2012 mm was a remarkable success of the Standard Model 
and established that the Higgs boson is indeed responsible for the breaking of the electroweak 
symmetry ms]. In the context of supersymmetry (for a review see [6]), the observed Higgs boson 
can be identified as the light CP-even state (see [MS]). Further, the fact that the observed 
Higgs boson mass is ^ 126 GeV implies the need for a large loop correction to its mass to raise 
its tree value (which lies below Mz) to the desired experimentally observed value. A large loop 
correction in turn implies a large supersymmetry-breaking weak scale on the order of a few TeV. 
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However, a high SUSY scale does not preclude some of the sparticle masses lying much lower in a 
widely split sparticle mass spectrum. In this work we focus on the possibility that, in spite of the 
overall SUSY scale being high, there might exist light squarks, specifically light stops, which may 
be discoverable at LHC RUN-IL 

In models with large trilinear couplings, a large split between the stop masses is automatic, and 
for the case when the mass gap between the light stop and the LSP is small, i.e., 

-m^o) « m^o (1) 

the decay of the stop is likely to produce soft jets, making detection more difficult. This possibility 
becomes a requirement if the model is constrained to satisfy WMAP US] and Planck m relic 
density constraints. Thus a satisfaction of the WMAP and Planck relic density constraints in this 
case requires that the lightest stop mass be within 20% or less of the LSP mass (for some recent 
works on light stops see [T8]-[28]). In the analysis of the light stops one must take account of the 
constraint that there be no instability arising from color and charge breaking minima. For most 
recent works on the necessary constraints for stability see [291130]. 


We will work within the framework of supergravity grand unification with radiative breaking 

of the electroweak symmetry (for a review see [35|). For this class of models it is predicted that the 
upper limit on the Higgs boson mass lies below 130 GeV (MSIESI. While the general framework of 
supergravity grand unification allows for non-universalites in the scalar boson sector, the gaugino 
sector, and the sector with trilinear scalar couplings, we will focus exclusively on the universal case, 
i.e. mSUGRA. The parameter space of SUGRA models has been analyzed in a number of previous 
works [36IH6] . Here we will constrain the parameter space so that the Higgs mass in the model is 
consistent with the LHC data and the relic density is consistent with the WMAP and Planck data. 
The resulting parameter space consistent with these constraints requires that the light stops lie is 
a narrow corridor between the LSP mass (which in all cases in this analysis is a neutralino) and 
1.2 times the LSP mass. 

In this region the relic density is satisfied by stop-neutralino coannihilatio gniiMB]. For this 
class of models, the mass difference between the stop and LSP is small and the decay 

products of the stop are relatively soft. It is found that the lightest stop can have a mass as low as 
375 GeV and stop pair production cross-sections can be as large as 2 pb. Further, using optimized 
signal regions it is found that the 375 GeV stop can have a ba discovery with 60 fb“^ of integrated 
luminosity. Additionally it is found that in this model stop masses as large as 600 GeV or larger 
with only a ^ 40 GeV mass gap between the stop and the neutralino can be discovered with the 
design parameters of the LHC RUN-H. 

The outline of the rest of the paper is as follows: In section we discuss the parameter space of 
interest, examining the stop-neutralino coannihilation region of mSUGRA parameter space under 
the constraints of the Higgs boson mass, WMAP and Planck relic density, and LHC RUN-I exclusion 
plots on the sparticle mass spectra. In section|^we carry out a signature analysis of a representative 
set of parameter points for the LHC RUN-H energy of 14 TeV. Here we analyze various decay 
channels and signal regions to determine the best avenues for the discovery of this class of models 
and determine the minimum luminosity needed for the discovery of each of the parameter points 
discussed. In section we investigate the direct detection of dark matter in models of the type 
discussed in section Conclusions are given in section 
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2 Parameter Space for Light Stops in mSUGRA 


An analysis of the parameter space of mSUGRA was carried out to identify regions with light 
stops that also satisfy the Higgs boson constraint, the relic density constraint, and the sparticle 
mass lower limits from LHC RUN-I m- The imposition of these constraints drastically reduces 
the parameter space of models in a manner discussed further below. The sparticle spectrum was 
generated using SoftSusy 3.5.1 [SHj with the mSUGRA input parameter set 

mo, Ao, mi, tan/3, sgn(/i) (2) 
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where mo is the universal scalar mass, mi /2 is the universal gaugino mass, and Aq is the universal 
trilinear scalar coupling (all at the grand unification scale); tan/3 = {H 2 )/{Hi)^ where H 2 gives 
mass to the up quarks and Hi gives mass to the down quarks and the leptons, and sgn(/i) is the 
sign of the Higgs mixing parameter which enters in the superpotential in the term iiHiH 2 > The 
analysis of the relic density was done using Micromegas 4.1.7 [59] and SLHA-format data files 
were processed using PySLHA [60] . 


2.1 Stop-Neutralino Coannihilation Region 


Satisfaction of the relic density constraint with a light stop requires stop-LSP coannihilation. The 
condition of coannihilation is typically 


Am 


m~o 

A1 




m^o 

A1 


< 0.2 


(3) 


Indeed, eq. ?? strongly hints that stop-neutralino coannihilation will be the dominant mechanism 
for satisfying the relic density constraint. It is thus of interest to map out the parameter space 
where the constraint of eq. ?? is satisfied. This is exhibited in the two panels of fig. [T] The left 
panel of fig. gives a three dimensional plot of the coannihilation region with mo, Ao/mo, and 
mi 12 as the three coordinates. The vertical axis on the right gives the color coding for Ao/mo- 
The analysis shows that the coannihilation region stretches out quite far in mo and mx /25 with 
mo extending past 20,000 GeV and mi /2 getting as large as 6000 GeV. The section of the curve 
where mo and mi /2 become very large is part of the hyperbolic branch of radiative breaking of the 
electroweak symmetry |6H - I64| . The right panel of fig. gives the projection of the left panel in the 
mo — mi 12 plane. The vertical axis on the right gives the color coding for the mass gap between 
the stop and the neutralino. The analysis of fig. was done under the Higgs boson mass constraint 
126 ± 2 GeV and the mass gap constraint defined in eq. 3, but no relic density constraints. 

It is possible to satisfy the WMAP and Planck relic density constraints in the coannihilation region 
shown in fig.[^ We consider two sets of relic density constraints: (i) the weak relic density constraint 
^LSP < 0.12, and (ii) the strong relic density constraint, which we take to be 0.0946 < f^LSP < 
0.1306, where the range is ±5 (j around the mean WMAP result. The weak relic density constraint 
allows for multi-component dark matter, i.e. that the relic density is made up of one or more dark 
matter components other than the neutralino, while the strong relic density constraint strictly 
requires that the dark matter is constituted only of the neutralino LSP. The range ±5 (j around the 
mean of WMAP is taken to allow for possible uncertainties in theoretical computations of the relic 
density. We give now an analysis of the parameter space for these cases. 
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Figure 1: Left panel: A three dimensional plot of the stop-neutralino coannihilation region satisfying 

eq. ?? and the Higgs boson mass constraint with three axes chosen as mo, Ao/mo, and mi and 

2 

with tan/3 = 10. The image is colored by Ao/mo values in the range [—2.2,—2.3]. Right panel: 

Two dimensional projection of the left panel in the mi — mo plane. The allowed regions are colored 
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by the stop-neutralino mass gap Am = m^ — m ~o. 


In the top two panels of fig. we investigate the parameter space for the case where we impose the 
weak relic density constraint (i). The axes in fig. |^are as in fig. Here one finds that the major 
difference between fig. and the top panels of fig. is that the coannihilation region has shrunk 
and the mass gap between the stop and the neutralino is reduced to < 45 GeV. In the bottom two 
panels of fig. we investigate the parameter space for the case where the stronger relic density 
constraint (ii) is imposed. The two-sided constraint in this case narrows the allowed range of the 
mass gap between the stop and the neutralino to between 30 and 40 GeV, as exhibited by the 
vertical bar to the right of the bottom right panel. As in fig. the analysis of all panels in fig. 
is done under the Higgs boson mass constraint. Similar plots arise for other tan/3 values, such as 
tan /3 = 30 and tan /3 = 50. 

The analysis of figs. and [fallows for a few conclusions. While fig. [^suggests that the coannihilation 
region as defined by eq. ?? continues in mo past ^ 25 TeV, there is an upper cutoff around 
mo ^ 8000 GeV once the relic density constraint is taken into account, as shown in fig.|^ Figure 
displays the ii — Xi mass plane and demonstrates that once the relic density constraint is applied, 
the allowed mass gap is greatly constrained to be between 30 and 40 GeV. From this analysis 
one concludes that once relic density is considered, the mass gap always lies below mt 

so that the on-shell decay ii tx\ does not occur, and the dominant decay for the stop in the 
coannihilation region is ii cxi- This decay remains dominant in the region where Am^^ has 
an upper limit of mw + m^ ~ 85 GeV. Of course the off-shell decay ii W^bxi can still occur, 
which will produce signatures of the type ii bff'xi^ where //' arises from the decay of the VF*. 
The jets arising from the decay will be soft, but high Pt jets could arise from initial and final state 
radiation (ISR and FSR). 
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Figure 2: Top panels: Same as the two panels of fig. [^except that relic density constraint f^LSP < 
0.12 is imposed. Bottom panels: Same as the two panels of fig. [^except that the strong relic density 
constraint 0.0946 < f^LSP < 0.1306 is imposed. 



Figure 3: Left panel: A plot of the neutralino mass versus stop mass for the full coannihilation 
region defined by eq. ?? and displayed in fig. The image is colored by the mass gap Am between 
the stop and neutralino masses. All masses are in GeV. Right panel: Same as left panel but 
restricted to only those points which satisfy the strict relic density constraint (ii) and which are 
discoverable with a luminosity oi C < 3000 fb“^, according to the analysis in section 
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3 Signal Analyses for Light Stop Models in LHC RUN-II 


To assess the appropriate channel for detecting a supersymmetric signal arising from the coannihi¬ 
lation region of parameter space allowed by the Higgs boson mass constraint and the WMAP and 
Planck relic density constraints, a Monte Carlo analysis was performed for a representative subset 
of coannihilation region points. The analysis was conducted using MadGraph 2.2.2 [65H67] . which 
includes Pythia 6.4.28 [6H] for hadronization of decays and Delphes 3.1.2 [69] for detector sim¬ 
ulation and event reconstruction. MadGraph calculates the possible decay diagrams for a requested 
set of decays. From there, the MadEvent submodule of MadGraph takes as inputs the sparticle mass 
spectrum, decay widths, and branching ratios, and uses this information to find the decay cross 
sections of the chosen decays. With these cross sections MadEvent can simulate decay events, using 
Pythia to complete the hadronization of resulting color-charged particles. Once event generation 
is complete, Delphes is used to reconstruct events within the geometry and triggering environment 
of a particular detector, outputting a file formatted for later analysis with ROOT 5.34.21 IZO!. 

For each simulated point of mSUGRA space, the sparticle spectrum was calculated with SoftSusy 
and input to SUSY-HIT 1 .5 m, which uses SDECAY |72l [73] and HDECAY [TU [75] to calculate the 
decay widths and branching ratios of the sparticles from the sparticle spectrum. Particular to the 
coannihilation region spectra studied here, the flavor-violating stop decay module of SUSY-HIT [76] 
is used in this analysis to obtain more accurate results for the compressed spectrum stop decays. 
The results calculated by SUSY-HIT in flavor-violating stop decay mode are inserted back into the 
general spectrum output given by the vanilla SUSY-HIT code. 

To ensure that every possible decay is considered, Feynman diagrams for the decays p p ^ SUSY 
SUSY were generated in MadGraph, where “SUSY” represents any MSSM particle. Both ISR and 
FSR are included in the analysis. These diagrams were then used to generate 50,000 MSSM decay 
events for a collection of points in the coannihilation region of mSUGRA space (see table [^, using 
each point’s calculated sparticle mass spectrum, decay widths, and branching ratios, as described 
above. For each mSUGRA parameter point, MadGraph generated an overall SUSY decay cross 
section, as well as cross sections for each individual decay mode, and then randomly generated 
decays according to those cross sections, hadronizing jets using Pythia. Once all 50,000 events were 
generated, Delphes completed the detector simulation using the ATLAS detector file, packaging 
the results into a ROOT output file. 

Search analyses were performed on the generated event set for each mSUGRA parameter point. 
The analyses were written as ROOT scripts implementing the search region details for previously 
published searches at 8 TeV. Specifically, the low-multiplicity jets plus missing transverse energy 
search signal regions from m and the monojet search signal regions from m were implemented 
and performed on the simulated MSSM decays. 

To compare to the background, the analyses were also performed on the pre-generated SNOWMASS 
published backgrounds m- From the total MSSM cross section calculated by MadEvent and 
the number of events simulated, an implied luminosity for each point was calculated to allow 
direct comparison to the backgrounds, which were also scaled individually by implied luminosity. 
The many background processes which comprise the SNOWMASS background set are scaled and 
combined in fig. Each individual process is represented by a different color in these plots. The 
various background samples are grouped according to the generated final state, with a collective 
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Figure 4: Top panels: Full SNOWMASS SM background after triggering cuts and a pre-cut 
of > 100 GeV, broken into individual processes and scaled to 100 fb“^. Left panel gives 


MEfr(incL) and right panel gives Individual data sets are labeled according to eq. (|^. Bot¬ 
tom panels: Meaf( incl.) (left) and (right) for the full SNOWMASS SM background, scaled to 

100 fb“^, after cuts which define the 2jl-opt signal region (see table and section 3.2) are applied. 


notation given by 

J = {i/, (i, J, 5, s, c, c, 6,5} , 

L = {e+,e,/x+,/^,r+,r,i/e,t'/„^'T} , 

1/ = {t^+,t^,Z,7} , (4) 

T = {t,i} , 

H = {ho} . 

in general, events with gauge bosons and SM Higgs bosons in the final state are grouped into a 
single “boson” (B) category. Thus, for example, the data set “Bjj-vbf” represents production via 
vector boson fusion of a gauge boson or a Higgs with at least two additional light-quark jets. 

Signal region cuts are applied to these backgrounds in addition to SUSY Model signals, and then 
the backgrounds for each process are summed before calculation of the integrated luminosity for 
discovery. The top panels of fig. [^illustrate two key kinematic quantities, Meaf and F™^®® for the 
full background dataset, after minimal cuts for trigger simulation and a F™^®® > 100 GeV pre-cut. 
Here and throughout, Meaf is specifically MEAF(incl.), defined as the scalar sum of F^®® and the 
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Model 

mo 

Ao 

mi 

2 


Al 


rrih 


m-g 

m^o m~± 

A2 Xi 

i. 

2910 

-6547 

770 

374 

337 

37 

124.5 

2316 

1865 

652 

ii. 

2990 

-6727 

790 

383 

346 

37 

125.6 

2376 

1900 

669 

hi. 

3150 

-7087 

830 

403 

365 

38 

127.9 

2497 

1989 

704 

iv. 

3350 

-7570 

900 

436 

397 

39 

124.8 

2663 

2142 

764 

V. 

3470 

-7842 

930 

447 

410 

37 

127.1 

2754 

2208 

790 

vi. 

3730 

-8467 

1020 

490 

452 

38 

124.1 

2968 

2403 

868 

vii. 

3810 

-8648 

1040 

498 

461 

37 

125.3 

3029 

2447 

885 

viii. 

3920 

-8898 

1070 

515 

475 

40 

125.5 

3114 

2512 

911 

ix. 

4150 

-9420 

1130 

543 

503 

40 

127.8 

3289 

2644 

964 

X. 

4350 

-9918 

1210 

579 

540 

39 

124.1 

3462 

2813 

1032 

xi. 

4460 

-10168 

1240 

595 

554 

41 

124.4 

3547 

2878 

1059 

xii. 

4610 

-10510 

1280 

614 

572 

42 

125.4 

3692 

2965 

1093 


Table 1: The subset of mSUGRA parameter points which satisfy the Higgs boson mass constraint 
and the relic density constraint in the stop coannihilation region and which were chosen for LHC 
analysis. The mSUGRA parameters are given (all points have tan/3 = 10 and /i > 0), followed by 
key superpartner masses (in GeV). Our simulations show that the stop-neutralino mass combina¬ 
tions arising from Models i-xii lie outside of the exclusion plots of RUN-I of the LHC. 


Pt of all jets with prU) ^ 40 GeV. The lower panels demonstrate the same quantities in the full 
background after applying the specific cuts which define signal region 2jl-opt, defined later in the 
text (see section 3.2). 


3.1 LHC Production and Signal Definitions 

In table we present a set of model points in the mSUGRA parameter space with light stops which 
give a Higgs boson mass consistent with LHC measurements and have a relic density consistent 
with WMAP and Planck observations. For the sake of consistency, all points in table [T] assume 
tan/3 = 10, though the results we will discuss are largely insensitive to the precise value of this 
parameter. Our simulation of these points at ^/s = 8 TeV shows that they would have escaped 
detection using the total integrated luminosity accumulated at LHC RUN-I [57]. In all cases, for 
the parameter points listed in table the light stop is the NLSP and its dominant decay mode is to 
a light chargino and a charm quark, which is a flavor-changing process. The other lightest particles 
are the second lightest neutralino the light chargino and the gluino g. Thus in the model 
points of table the sparticle mass hierarchy i^ 

m~o < mi^ < m~o < m~± < rrig . (5) 

Each of the NNLSP decays, i.e. the decays of X2^Xi’> involves a light stop in their dominant 

decay. However, the decay of ti in the coannihilation region will yield soft jets because of the small 
mass gap between its mass and the LSP mass, which makes it difficult to observe it directly. 

In tablewe give the production processes involving U, X 2 ^ 9 RUN-H of the LHC for 

each of the parameter points listed in table The dominant production modes generally consist 


^The hierarchy eq. ?? is mSP[tla] in the notation of Table 5 of [8^ . 








Model 

gg^^ 

iiii 

V 

hh 

qq 

xlxi 

qq 

xtxi 

gg gg 

qq gg 

BR(txtx — V 

cx\) 

i. 

2.0 

0.25 

1.0x10-2 

4.8x10-3 

9.0x10-^ 

2.1x10-^ 

89.7% 

ii. 

1.7 

0.22 

8.8x10-^ 

4.3x10-3 

6.9x10-^ 

1.8x10-^ 

88.0% 

hi. 

1.3 

0.18 

6.9x10-^ 

3.3x10-3 

4.1x10-^ 

1.2x10-^ 

83.9% 

iv. 

0.85 

0.12 

4.6x10-^ 

2.1x10-3 

1.8x10-^ 

5.9x10-3 

81.4% 

V. 

0.74 

0.11 

3.9x10-^ 

1.8x10-3 

1.3x10-^ 

4.4x10-3 

86.0% 

vi. 

0.44 

7.1x10“2 

2.4x10-3 

1.1x10-3 

3.8x10-3 

1.8x10-3 

80.9% 

vii. 

0.40 

6.6x10-2 

2.4x10-3 

1.0x10-3 

3.2x10-3 

1.5x10-3 

83.4% 

viii. 

0.33 

5.6x10-2 

1.9x10-3 

8.9x10-^ 

2.1x10-3 

1.1x10-3 

74.7% 

ix. 

0.24 

4.3x10-2 

1.4x10-3 

6.4x10-^ 

9.4x10-3 

6.3x10-3 

71.0% 

X. 

0.17 

3.1x10-2 

9.1x10-^ 

4.3x10-^ 

3.9x10-3 

3.0x10-3 

71.9% 

xi. 

0.14 

2.7x10-2 

7.9x10-^ 

3.7x10-^ 

2.9x10-3 

2.2x10-3 

61.5% 

xii. 

0.12 

2.3x10-2 

6.6x10-^ 

3.0x10-^ 

1.6x10-3 

1.5x10-3 

60.9% 


Table 2: Production cross sections in pb for LHC RUN-II of the dominant supersymmetric modes 
gg hh, qq ^ iiiu qq xlxi^ qq xixi, gg gg, and qq gg, as well as the branching 
ratio for the process iiii —> cx^, for the set of parameter points given in table 


of the following set: 


gg iiii, qq hh, qq X2X1 , qq^xtxi, gg^ M, and qq^M- ( 6 ) 

The largest production cross section is for the stops, followed by the production of the weak 
gauginos, and then the gluino. For the parameter regime presented in table the typical mass 
difference ~o is consistently near 40 GeV, with the ffavor-violating decay ii cxi typically 

representing 60% or more of the total decay width of the lightest stop. 

Given the dominance of stop production followed by ii our signature analysis concentrates 

on two ATLAS searches: the low-multiplicity jets plus missing transverse energy search m 

and the dedicated stop search of Egg The multi-jet search defines fourteen signal regions with 
jet multiplicities between 2 and 6 jets, and with varying requirements on the inclusive effective 
mass MEff(inc.), designated as “loose” (low MEff(inc.)) to “tight” (high MEff(inc.)). The dominant 
signal for the points in table involves typically 2-4 reconstructed jets with prU) > 20 GeV and 
a low effective mass. Examination of all 14 signatures therefore reveals that the two-jet “loose” 
signature (2jl) is the most effective topology for searching for these models. This signal region has 
the selection requirements as listed in the left panel of table 

The dedicated stop search of m involves two types of topology. The first is a monojet-like sig¬ 
nature, divided into three signal regions—Ml, M2, and M3—defined by increasingly stringent re¬ 
quirements on the PtUi) of tho leading jet and the missing transverse energy. The second topology 
also involves a hard leading jet, with priji) requirements similar to that of the monojet analysis, 
but also requires a charm tag on at least one of the sub-leading jets. The multi-variate technique 
that generates the charm tag is not something that is easily reproduced, nor modeled in Delphes. 
We therefore focus only on the mono jet searches in this work. The most successful mono jet search 
for this parameter region is Ml, which has the selection requirements specified in the right panel 
of table m 

^An alternative set of signal regions for supersymmetry discovery are used by the CMS Collaboration, see e.g., 
[ 801 IST] which involves razor variables. 
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Requirement (2jl SR) 

Value 



(GeV) 

>160 

Requirement (Ml SR) 

Value 

PtUi) (GeV) 

>130 

(GeV) 

> 220 

PtU2) (GeV) 

> 60 

PtUi) (GeV) 

> 280 

A<^(jeti,2,^r")mm 

> 0.4 

AcPiiet,Ef^^U^ 

> 0.4 

(GeV^/^) 

> 8 

At most 2 other jets with pr (GeV) 

> 30 

MEff(inc.) (GeV) 

> 800 




Table 3: Left: The selection criteria used for the signal region 2jl, where “2j” stands for two jets 
and “1” stands for loose in the nomenclature of Table 2 of the ATLAS analysis This signal 
region is found to be the most dominant for the parameter set given in table [Tj Right: The selection 
criteria for the monojet signal region Ml from the ATLAS dedicated stop search m- This signal 
region is subdominant for the parameter set of interest. 


Model 

jC for 5a 
discovery in 
2jl-opt SR 

jC for 5a 
discovery 
in 2jl SR 

jC for 5cr 
discovery 
in Ml SR 

i. 

61 

88 

369 

ii. 

79 

112 

396 

hi. 

100 

135 

535 

iv. 

148 

202 

739 

V. 

180 

242 

111 

vi. 

322 

446 

1561 

vii. 

367 

522 

1881 

viii. 

502 

706 

2460 

ix. 

725 

1025 

4035 

X. 

1214 

1716 

6118 

xi. 

1621 

2248 

8539 

xii. 

2140 

2994 

10538 


Table 4: Analysis of the discovery potential for supersymmetry for the parameter space of table [T| 
using the leading channel 2jl, the sub leading signal channels Ml, and the optimized channel 2jl-opt, 
where the minimum integrated luminosity needed for 5a discovery is given in fb“^. 
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Mg„ for Model i in M1 at 369 fb ' e 7“ for Model i in Ml at 369 fb ' 
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Figure 5: Left panel: Distribution in Meaf for the Ml signal region for Model i. Plotted is the 
number of counts for the SUSY signal and the square root of the total SM backgrounds. The 
analysis is done at 369 fb“^ of integrated luminosity, which gives a ba discovery in this signal 
region. Right panel: The same analysis as in the left panel but for 


Using the techniques described above, we analyzed each of the model points of table to identify 
a minimum required luminosity for ba signal discovery in the signal regions of interest. In 

^ ^ Vbackground ^ o o 

the analysis here, only statistical errors were considered and no attempt was made to estimate 
systematic errors. Such an effort requires detailed and specialized knowledge of the experimental 
apparatus. Furthermore, systematic errors at ^/s = 14 TeV are expected to be very different in 
form and magnitude than the published systematic errors at ^/s = 8 TeV and at lower luminosity. 
Our estimates of the integrated luminosities given here are based on single-channel analyses, and 
one expects that combining channels can potentially result in discovery with even smaller integrated 
luminosities, even after accounting for systematic error. 


In addition, because the tightly constrained parameter space yields a similarly constrained sparticle 
spectrum and signal, variations of these signal regions were tested in order to optimize them for 
sparticle spectra with the particular properties of tables and (see section 3.2). The minimum 
integrated luminosity results of these analyses are presented in table It is found that these 
model points will be discoverable with integrated luminosities beginning at ^90 fb“^ using existing 
searches or ^60 fb“^ using an optimized search, described in section [3^ 


Results for the MI signal region are represented in fig. which displays the Model i signal and the 
square root of the summed background after MI cuts. The left panel plots Meaf while the right 
panel plots These figures are shown at 369 fb“^, the necessary integrated luminosity for ba 

discovery for Model i using signal region MI. 


Likewise, a representative result for the 2jl signal region is shown in fig. Like fig. this displays 
the signal for Model i superimposed upon the square root of the total summed background, both 
after cuts, this time in 2jl. The left panel plots Meaf and the right panel plots The integrated 

luminosity used is 88 fb“^, which gives ba discovery of Model i using signal region 2jl. 


Comparing the results of table for the dominant signal and for the subdominant signals, we find 
that the integrated luminosity needed for the subdominant signal MI can be larger by a factor 
of 3 or more. All other signal regions from and m that were analyzed required even higher 
integrated luminosities for ba discovery. This verifies the claim that the 2jl is the dominant signal 
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Figure 6: Left panel: Distribution in Meap for the 2jl signal region for Model i. Plotted is the 
number of counts for the SUSY signal and the square root of the total SM backgrounds. The 
analysis is done at 88 fb“^ of integrated luminosity, which gives a ba discovery in this signal region. 
Right panel: The same analysis as in the left panel but for 


for discovery of the models listed in table However, because 2jl is a general-purpose search 
optimized at 8 TeV, it is possible to improve upon its performance for this parameter region. 


3.2 Optimizing the Signal Regions 

As mentioned above, because the leading signal region 2jl is a general-purpose signal region de¬ 
veloped for application at 8 TeV, it should be possible to improve on its performance by closely 
examining the cuts made to certain key parameters. Three primary changes were identified for 
analysis and comparison to the baseline signal regions: Relaxing the kinematic cuts on MEff(inc.), 
or PtU)^ increasing the cut on Ac/) between and the two leading jets, and making an 

additional cut on a new parameter given by the ratio pT{ji)/E^'^^^. Optimizations on these pa¬ 
rameters were performed on the 2jl and Ml signal regions described in table as well as on signal 
region 3j from izzi, which differs from 2jl primarily by requiring a third jet with prijs) > 60 GeV, 
and by imposing a much tighter cut Meaf > 2200 GeV (as opposed to Meaf > 800 GeV). 

With values of ~o of around 40 GeV, there is only energy available in the ii cxi to make 

a fairly soft jet—ISR is generally relied upon to produce high-p^ jets in the final state. Thus, it 
is reasonable that relaxing the kinematic cuts on MEAF(inc.) could yield an improved signal region. 
For the case of 2jl, reducing the MEAF(incl.) cut from 800 GeV to 400 GeV improved the result 
from the base case by almost a factor of 2 for the lightest points in the parameter space. But 
as parameter points become heavier this advantage decreases, and eventually the modified signal 
region becomes worse. This trend is reflected also if the kinematic cuts in Ml are reduced so 
that both > 150 GeV and priji) > 150 GeV are required; performance of the signal region 

improves for the lightest points but becomes worse for higher-mass regions of the parameter space. 
Reducing the Meaf (inch) cut on 3j from 2200 GeV to 800 GeV, a value that matches the cut for 
2jl, dramatically improves the 3j performance, leaving it only slightly worse than 2jl. In fact, after 
this change there is little that distinguishes the two signal regions; the modified 3j requires a third 
jet with a minimum prijs) > 60 GeV, while 2jl does not look for a third jet. Thus, 3j is dropped 
from consideration because its optimization evolves it into 2jl. 


12 








A4,(j ^, £™“) at 100 fb’’ 4^02, at 100 fb'^ 



0 0.5 1 1.5 2 2.5 3 0 0.5 1 1.5 2 2.5 3 . 


Figure 7: Left panel: Distribution of A(/)(ji, against the square root of the full SM background 

at 100 fb“^ for Model i before any cuts are applied. Right panel: The same as in the left panel but 
for 

Because the signal is expected to consist primarily of a single jet recoiling against a higher 

A(/)(j, is expected in the signal relative to the more uniform background. This is demon¬ 
strated in fig. where the left panel shows the shape of A(/)(ji, against the square root of 

the summed background for Model i before any cuts are applied, and the right panel shows the 
same for A(/)(j 2 , Based upon these results, Ac/) cuts were increased from 0.4 to | for ji and 

from 0.4 to 1.0 for j 2 . Making this change improved the 2jl results by ^30%, yielding the signal 
region described in table as 2jl-opt. Making this same change in the Ml signal region yielded a 
poorer search; the integrated luminosity required for 5a discovery increased. 

Finally, motivated by recent suggestions in [8311H1] ,we considered including an additional kinematic 
parameter, r = pT{ji)/. Because this parameter measures the degree to which the recoiling 
missing energy is concentrated in a single jet, it was expected that it would help distinguish this 
signal. However, when r > 0.5 was applied to 2jl and Ml, with and without optimization in Ac/), 
the new signal region was worse in every case. 

After these investigations, the final optimal signal region was deemed to be 2jl-opt, which, as 
described above, is the same as 2jl (see table but with Ac/) cuts increased from 0.4 to | for ji 
and from 0.4 to 1.0 for j 2 . This signal region boasts a consistent 30% increase in performance as 
compared to 2jl across the full range of the parameter space. 

In fig. we exhibit Meap and distributions in the optimized signal region 2jl-opt for three 

selected models (i, vii, and xii). Plotted are the number of counts per 30 GeV energy bins for 
the SUSY signal and the square root of the full Standard Model backgrounds. In the left panel 
of each figure we give the analysis for Meap and in the right we give the analysis for The 

top panels show the analysis for Model i. Here the distribution is given at 61 fb“^ of integrated 
luminosity, which is found to be minimum integrated luminosity which gives a 5cj discovery signal 
for supersymmetry in 2jl-opt. Identical analyses are given in the middle panels for Model vii, but 
at 367 fb“^ of integrated luminosity, which gives a 5cj discovery signal for that point, and in the 
bottom panels for Model xii. Here one finds that a 5cj discovery requires a minimum of 2140 fb“^ 
of integrated luminosity. Thus Models i, vii, and xii are all eventually discoverable with the design 
parameters of the machine at y/s = 14TeV. 


13 











Mg„ for Model i in 2jl-opt at 61 fb'^ 


(O 



800 900 1000 1100 1200 1300 1400 


Me„ (GeV) 

Mgjj for Model vii in 2jl-opt at 367 fb'^ 



800 900 1000 1100 1200 1300 1400 

Me„ (GeV) 

Mg„ for Model xii in 2jl-opt at 2140 fb'^ 


£ 



800 900 1000 1100 1200 1300 1400 


Me„ (GeV) 


for Model i in 2jl-opt at 61 fb‘^ 





0 200 400 600 800 1000 1200 1400 

E™'“ (GeV) 

Ej'^ for Model vii in 2jl-opt at 367 fb‘^ 



E7'“ (GeV) 


for Model xii in 2jl-opt at 2140 fb'^ 


£ 



0 200 400 600 800 1000 1200 1400 

E?*'** (GeV) 


Figure 8: Distribution in Meap (left panels) and (right panels) for the optimized two-jet 

signal region 2jl-opt. Top panels show Model i normalized to 61 fb“^ of integrated luminosity, 
while the middle and bottom panels give Models vii and xii, normalized to 367 fb“^ and 2140 fb“^ 
of integrated luminosity, respectively. These numbers give a ba discovery in this signal region for 
the three model points. Plotted in all three cases is the number of counts for the SUSY signal and 
the square root of the total SM background. 
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4 LHC and Dark Matter 


Interesting connections between LHC physics and dark matter have been noted in previous works 
(see [85H9T]). In the models considered here this connection is even stronger because the requirement 
that dark matter satisfy the WMAP and Planck relic density constraints forces the stop and 
neutralino masses to lie within ^ 30 GeV of each other for the lightest stops, and restricts the 
allowed parameter space of the model in the mo — mi /2 plane to a very narrow strip. Further, since 
^o/'^o is also determined to be in a very narrow range centered at ^ —2.25 in order to generate 
the right Higgs boson mass, the model is very predictive. It is of interest then to investigate 
the implications of these constraints for the direct detection of dark matter. For the parameter 
space of the model, the neutralino turns out to be almost exclusively a bino with extremely small 
wino and Higgsino content. Therefore proton-neutralino scattering cross sections are expected to 
be small. The analysis of spin-independent and spin-dependent proton-neutralino cross sections 
is presented in table The spin-independent cross sections lie in the range between 10“^^ and 
10“^^ cm“^. While cross sections of this size are decidedly small, they could still be visible in 
the next generation LUX-ZEPLIN (LZ) dark matter experiment, which is projected to reach a 
sensitivity of 10 ^ Regarding the spin-dependent proton-neutralino cross section, 

the LUX-ZEPLIN will have a maximum sensitivity of 10“^^cm“^, which is still about three orders 
of magnitude smaller in sensitivity than what is needed to observe the spin-dependent proton- 
neutralino cross section exhibited in table Thus the observation spin-dependent cross sections 
will be more difficult. 
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Model 

pxV 

X 10^® 
pxV 

i. 

3.5 

9.6 

ii. 

3.2 

8.6 

hi. 

2.6 

6.9 

iv. 

2.6 

5.2 

V. 

2.2 

4.5 

vi. 

2.1 

3.2 

vii. 

1.9 

3.0 

viii. 

1.8 

2.6 

ix. 

1.5 

2.1 

X. 

1.6 

1.7 

xi. 

1.5 

1.5 

xii. 

1.4 

1.3 


Table 5: CDM proton-neutralino spin-independent and spin-dependent cross sections 

in units of cni“^ for the set of model points of table [l] None of the parameter points are ruled out 
by the current dark matter experiments, in particular by LUX. 

5 Conclusion 

In this work we have analyzed the lightest stop masses that can arise in high-scale models consistent 
with the Higgs boson mass constraint, relic density constraints from WMAP and Planck, and 
constraints on sparticle mass spectrum from RUN-I of the LHC. The specific focus of this work was 
mSUGRA, where it is found that stop masses as low as 400 GeV or lower can still exist consistent 
with these constraints. However, in the model space analyzed it is found that the stop must be the 
NLSP, with a mass ^1.1 times the LSP mass. This restriction is needed in order to satisfy the 
relic density constraint via stop-neutralino coannihilation. 

Requiring that the lightest neutralino supplies all of the thermal relic abundance then puts a lower 
bound on Am = — m ~o of approximately 30 GeV. Over the range of parameter space accessible 

in the lifetime of the LHC (i.e. within 3000 fb“^ of integrated luminosity), this gap never exceeds 
45 GeV. Thus neither the decay ii txi nor on-shell VP-decays of the form ii Wbxi occur. 
Rather the dominant decay of the stop in the allowed parameter space is the flavor-violating process 
ii cxi- 

Production of superpartners is dominated by stop pair production, with all other processes reduced 
by approximately two orders of magnitude. The analysis of this work shows that the discovery of 
supersymetry for the class of models with light stops discussed here can occur with an integrated 
luminosity as low as ^ 60 fb“^ at the LHC RUN-H, and will be dominated by searches involving low 
jet multiplicity, lepton vetoes, and large missing energy. We have carried out a signature analysis 
for this class of models and find the dominant signature to be 2jl, where the characteristics of the 
2jl SR are given in table We have also carried out an analysis of the subdominant signature Ml 
and an optimized signature 2jl-opt, the characteristics of which are given in table It is shown 
that the light stop can be discovered in 2jl with 88 fb“^ and in 2jl-opt with 60 fb ^ of integrated 
luminosity at the LHC RUN-H, while a confirmation of the subdominant signal will require three 
times as much luminosity. The minimum integrated luminosity can no doubt be reduced with 
further improvement in the efficiency to tag charm-jets [94II96] . 
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We also discussed the dark matter associated with this class of models. It is shown that the spin- 
independent proton-neutralino cross section could be within reach of the next generation LUX- 
ZEPLIN dark matter detector while the spin-dependent proton-neutralino cross section will be 
more difficult to observe. Any observed signal will be comparable to, but slightly larger than, 
coherent scattering by atmospheric neutrinos m- 
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